Random Walks through the Ensemble: Linking Spectral Statistics with Wavefunction 

Correlations in Disordered Metals 
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We use a random walk in the ensemble of impurity configurations to generate a Brownian motion 
model for energy levels in disordered conductors. Treating arc-length along the random walk as fic- 
titous time, the resulting Langevin equation relates spectral statistics to eigenfunction correlations. 
Solving this equation at energy scales large compared with the mean level spacing, we obtain the 
spectral form factor, and its parametric dependence. 
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Statistical properties of the spectra of finite quantum 
systems have been a focus for research in three succes- 
sive contexts: nuclear physics 0], the semiclassical limit 
of quantum mechanics Q , and studies of mesoscopic con- 
ductors ||||]. A unifying idea is that the energies of in- 
dividual eigenstates are frequently neither calculable nor 
interesting: instead, the concern should be with eigen- 
value correlations, which typically are independent of 
many details of the Hamiltonian. 

Disordered mesoscopic conductors bring to this field 
both new behaviour and a natural ensemble for a statisti- 
cal description - the ensemble of impurity configurations. 
For weak disorder, new behaviour arises because there 
can be a broad window in time, and hence a correspond- 
ing energy interval, between the scale, t e \, for electron 
scattering from impurities and that for diffusion across 
the system, t CIg ~ L 2 /D (where L and D are the system 
size and diffusion constant). And at the mobility edge, 
specific, critical spectral statistics are expected [[3]||. The 
ensemble average provides the departure point for estab- 
lished perturbative Q and non-perturbative j|] calcula- 
tions of spectral correlations in disordered metals. 

In this paper we describe an alternative approach, in 
which the energy level distribution is averaged along a 
random walk through the ensemble. There are several 
precedents for study of levels as a function of position 
in the space of Hamiltonians. Most recently, a number 
of authors tffl, in particular Szafer, Simons and Altshuler 
pL have investigated parametric statistics: level correla- 
tions between different points on a smooth path in this 
space. Earlier, in connection with the semiclassical limit, 
Pechukas Pj used motion along such a path, with coor- 
dinate A, to generate a dynamics for the one-dimensional 
gas formed by levels on the energy axis. And originally, 
in the context of random matrix theory (RMT), Dyson 
JTfj| , employed a random walk through the matrix ensem- 
ble as the foundation for Brownian dynamics of levels, 
with arc-length, r, being the fictitious time. Pechukas' 
and Dyson's ideas are linked, as shown schematically in 
Fig. 1, by the usual relation JnJ between the end-to-end 
distance and the length of a random walk, A 2 ~ r. 



FIG. 1. Smooth [-ff(A)] and Brownian [H(t)] paths through 
the space of Hamiltonians. 

In previous work on level dynamics |p|-|l2|, eigenfunc- 
tion correlations have played no role, either (in the semi- 
classical limit) by assumption, or (for random matrices) 
by construction, RMT having no preferred basis. In con- 
trast, for disordered metals, the basis of position states is 
singled out: it is in this basis that the impurity potential 
is a diagonal operator. Further, specific eigenfunction 
correlations must be present if, for example, wavepackets 
spread diffusively. We show here that these correlations 
result in a novel Brownian level dynamics, from which 
we obtain the features of spectral correlations particular 
to mesoscopic conductors. 

Our results are expressed in terms of the quantum re- 
turn probability, p(t), which is required as input in this 
approach. Consider the time-evolution of a wavepacket 
that initially occupies a small volume Iq |13| ]: p(t) is 
the probability density to remain in this volume at the 
time t. For weak disorder, the wavepacket spreads diffu- 
sively, and p(t) is known from semiclassical arguments; 
at the metal-insulator transition, scaling ideas relate p(t) 
to fractal properties of eigenstates. 

We shall be interested in the two-level correlation func- 
tion (TLCF) and the spectral form factor. Let E n (r) be 
the energy levels of a one-parameter family of Hamiltoni- 
ans, H{t). The density of states per unit volume (DoS) 
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is p{E, t) = L d J2 n S(E — £ n (r)) , and the mean level 
spacing is A=(pL d )- 1 , with p={p(E,r)). The TLCF is 



R(E,r)=p- 2 (p(E,r)p(0,0))-l 



and the spectral form factor is 



1 r°° 

^(i,r) = — y e~ iEt ^R(E, r) dE . 



(1) 



(2) 



We obtain K(t, r) for times shorter than the Heisen- 
berg time, t H = Ti/A, in terms of p(t). We de- 
fer discussion of the parametric dependence until the 
end of this paper. The non-parametric form factor, 
K(t) = K(t,r= 0), is 



K(t) = 



1*1 K*) 



l+(TThp)-lf$ P (t>)dt> 



In the diffusive regime, t e \ < t < i e rgj at leading order 
reduces to the classical return probability for random 
walks, multiplied by a symmetry factor 2/(3, where = 
1 . 2 or 4 is the usual index (l41 : 



Po(t) 



(3(ATrDt) d / 2 ' 



(4) 



There exist quantum (weak-localization) corrections to 
Eq. (0), which are smaller by a factor g^ 1 , where go 3> 1 
is the dimensionless conductance at scale of the elastic 
mean free path, I. Noting that in the ballistic regime, 
t < t e \, pit) saturates at po(t c \) ~ one sees that the 
second term in the denominator of Eq. (^|) also gives weak 
localization corrections: it is of order (t c \A/h)(L/£) d ~ 
g^ 1 for d > 2, and of order g^ 1 \n(t/t c \) for d = 2. Ne- 
glecting these corrections in both the numerator and de- 
nominator of Eq. , it reduces to 



K (t) 



A 

2nh 



L d \t\ Po (t). 



(5) 



To leading order, this expression is also valid in the er- 
godic regime, t crg where the classical return 

probability saturates at (2/ (3)L d so that the second term 
in the denominator in Eq. (g) is of order t/t H <C 1 and 
may be neglected. 

Equation (|^) coincides with the result obtained in 
Rcf. [|l5|| , using the diagonal approximation in semiclas- 
sical periodic-orbit theory. The Fourier transform of 
this expression corresponds to the TLCF obtained orig- 
inally by Altshuler and Shklovskii from the lowest or- 
der of the diagrammatic expansion Q. In the diffusive 
regime, i?(s,0) ~ A c is d ^ 2 ^ 2 , and in the ergodic regime, 
R(s, 0) ~ — l/s 2 |l6| ], where s = E/A, and Ad is a numer- 
ical coefficient which is zero for d=2 Jlt| . 

There are indications that our approach is useful be- 
yond the diagonal approximation. We have checked Eq. 



([|) in the diffusive regime by direct diagrammatic expan- 
sion of both sides in powers of g^" 1 , calculating the first 
two non-trivial orders for d = 2. This is the most inter- 
esting case because, for d = 2, po(t) oc 1/t, so that K (t) 
is constant and R(s, 0) = in the diagonal approxima- 
tion. We find |ll| that Eq. (^) is exact up to order g^ 2 . 
For /3=1, 4 this means that both the leading (gQ 1 ) con- 
tribution to R(s, 0) and the g^~ 2 correction are exact. 
In the unitary case ((3 — 2) the first order terms in the 
numerator and denominator of Eq. (j^) cancel, in agree- 
ment with the known absence of corrections to Ko(t) at 
this order, and it is the g$ 2 contribution which governs 
R(s, 0); again it is given exactly by Eq. (||). For d > 2, 
first order corrections to the diagonal approximation are 
small by an additional factor, (i e i/i) d / 2_1 ; Eq. (||) repro- 
duces their form but not the numerical coefficient. On 
the other hand, at the mobility edge in d > 2, Eq. (||) is 
consistent with the results of independent calculations, 
as we discuss at the end of the paper. 

We next turn to the derivation of our results. Let the 
fictitious time r parameterize a Brownian path through 
an ensemble of disordered conductors, so that 



H{t) = H + f dr'y(r',r), 
Jo 



H = --V 2 + U(r) 



(6a) 
(6b) 



We take U(r) and V(r,r) to be Gaussian distributed 
with zero average and 



(U(r)U(r')) = (h/2Trpt el )6(r-r') , 
(V(t, r)V( T ', r')) = v 2 L d 5{t - r') 8{r-r') 



(7) 



We use perturbation theory to second order to calcu- 
late the change of E n (r) in response to V(T,r)ST: 



SE n = V nn + 



rmn*nrn 



Yn 



dr' (n\V(r',r)\m) , 



where (r\n) = ip n (T' ' ,r) are the corresponding eigen- 
functions of H(t'). We average over V noting that 



(V nm V mn ) = (V nn V m m) = v 2 J dr c„ m (r) , (8a) 
c nm (r)=L d fd d r\Mr,r)\ 2 \^ m (T,r)\ 2 . (8b) 



The level dynamics thus obeys the Langevin equation 



dE n (r) 2Sr c ra , n+ ;(r) 

IjtO 



dr 



E n {r) ~ E n+ i{ T ) 



where £ is a random force with (£) = and 



2 



(Ur)U(r')) =v 2 5(T-T')c nm (T). 



(9b) 



Both the spectrum of the random force and the drift term 
are expressed in terms of c nm (r), the eigenstate correla- 
tion function, which we later show to be proportional to 
the return probability p(t). 

The appearance of the eigenstate correlation function 
in Eqs. (||) is the essential difference between our Brown- 
ian level dynamics for disordered conductors, and that of 
Dyson for random matrices. The analogues to Eqs. 
in Dyson's work ||ic|| have the basis-independent form 



(10) 



In contrast to Eqs. (||), these display two simplifying fea- 
tures: they are independent of the eigenvectors, and they 
lead to diagonal random force correlations. As a result, 
the limiting level distribution in RMT can be calculated 
exactly, proceeding via a Fokker-Planck formulation [jl0| . 
For disordered conductors, the Langevin equation for en- 
ergy levels is not closed, since Eq. ( pa] ) involves c rml (r), 
and the random forces have off-diagonal correlations. 
In consequence, approximations are required, which are 
most transparent within the Langevin description. 

Two central assumptions are necessary in order to 
make progress. We believe that they are reasonable pro- 
vided one considers only behaviour at energy scales much 
larger than the mean level spacing. First, we replace 
Cnm{T~) by its average over the ensemble of Hq, which, 
within the energy window of interest, is a function only 
of the energy difference, E = E„ — E m : 



(c„ >n+ i(r)) = c{E) . 



(11) 



Thus we neglect correlations between fluctuations in 
Cnm{T) and those in p{E) pSj ], Second, we linearize Eqs. 
(^J) in the deviation of level separation E n — E m from its 
mean value (n — m)A, 

E n (r) = nA + e n (r). 

Then Eq. ( |9a| ) becomes 

tfe n (r) 



v 2 J2(e n+ i-e n )f(lA)+Ur) (12) 



where f(E) = dtp(E)/dE with <p(E) = c(E)/E. This is 
diagonalized by Fourier transform, using as coordinates 



A 



(13a) 



which, for < t <C ft-/ A, and e n t/h are essentially Fourier 
components of the DoS: 



£(t,r) 



2-Kpt 

In these variables, Eqs. (ft) become 



p(E, T)e- lEt/h dE . (13b) 



d£(t, t) 2-irhv 2 



F(t)£(t,T)+E(t,T) 



(14a) 



dr A 
<H(t, r) E(t', r')) = v 2 S(t + t')6(r-r')C(t) , (14b) 

where C(t), F(t) and E(t, r) are the Fourier transforms 
of c(E), f(E) and £„(r), defined as in Eqs. (|l|). For 
example, 



A 



iAlt/H 



c(E)e 



-iEt/h 



dE 
2nh 



It follows from the definition of f(E), after Eq. (|12|), that 
F(t) is related to C(t). Simple manipulations lead to 
F(t) = (t/Ti 2 ) J C(t')dt' , where we have used the fact 
that F(t) is an even function of t to fix the constant of 
integration. Finally, we note that, since c(E) = 1 for 
\E\ > h/t c \, C(t) contains a (^-function at t = 0, so that 



n*) = |Ui [ W c(t')dt>. 



(15) 



Now we solve Eqs. (|TJ), fixing at this point the units 
of t by setting v 2 — A/2irh: 

{ £ { t 1 r + r')£{t'y))=5(t+t')^^ ) e- F ^\ 

As £(t,r) is the Fourier transform of p(E,r), we imme- 
diately obtain the spectral form factor: 



C{t)A 
hj 4TrhF(t) ' 



-F(i)|T| 



(16) 



Note that the form factor is expressed in terms of eigen- 
function correlations via the interplay of the restoring 
force, F, and the noise correlator, C . It is the latter 
which is responsible for the difference between spectral 
statistics in the ergodic and diffusive regimes. 

We relate F and C to the return probability of a dif- 
fusing electron, by considering a wavepacket made from 
the eigenstates of H(t) and concentrated initially near 
the origin, in a volume of size £ d : 

a 

Here the summation is limited to J\f levels with energies 
E n fo/t c i so that N ~ fi/t c \A, and the normalization 
constant is A 2 = L'W -1 . The ensemble- averaged return 
probability, p(t) = (|*(0, t)\ 2 ), is given for t > by 

p(t) =A 2 J2 (c nm e^ E "' E - )t/h ) = 2nhpC(t), (17) 



where we have again (as in the solution of the linearized 
Langevin equation) neglected correlations between the 
fluctuations of {£u(t)} and c nm (r). 

Substituting for C(t) in Eq. [l6[ we find: 



3 



\t\p(t) 

M(t) 



exp 



M(t) 

2h 2 



tr 



i r 1 * 1 

M(t) = l + — p(t')dt' 



(18a) 
(18b) 



These expressions, which are the central result of this 
work, relate spectral statistics to the return probability 
and thus, via Eq. (|l7|), to the wavefunction correlations 
(|Sh|). For the non-parametric problem, Eq. (18a) reduces 
to Eq. (§), and the ensuing behaviour in the ergodic and 
diffusive regimes has been discussed following Eq. (Q). 

Parametric correlations in these regimes, from Eqs. 
(|l8|), are as follows. At leading order we ignore weak- 
localization contributions to pit) and M(t), approximat- 
ing them by po(t) and 1, respectively. Then, for t <C t H , 



Kit,T)=K (t)e-\ tT \/ 2h2 



(19) 



Its Fourier transform gives the parametric TLCF in both 
the ergodic regime (where it reproduces the envelope of 
the TLCF found in Refs. p|,^2[) and the diffusive regime 
(where it is consistent with the results of diagrammatic 
calculations). For E = we obtain 



i?(0, A) = 



A 4 , ergodic regime 

B c ig~ d / 2 \ d ~ i , diffusive regime 



(20) 



where A 2 = (ttt/TiA) and Bd is a constant which is non- 
zero in d = 2, so that the leading approximation is suf- 
ficient for the parametric TLCF in d = 2, in contrast to 
the non-parametric TLCF, R{E,T = 0). 

Finally, we turn to spectral correlations for a system at 
the metal-insulator transition |]||. At the critical point, 
power-law decay of the return probability is expected for 
t H , with the power related to a multi- 



id < t <C t e , 
fractal dimension of wavefunctions 20 



pit) oc t- 1+r '/ d . 



(21) 



where < r/ = d — c?2 < d, and di is the dimension asso- 
ciated with | -0| 4 pi] . In consequence, for the same range 
of t, Kit) is constant, a result obtained previously from 
a diagrammatic analysis ||. As in d = 2, the behaviour 
of RiE,r — 0) at the mobility edge is, in fact determined 
||, not by this limiting behaviour for K(t), but by the 
way in which the limit is approached as t/t H — + 0; we de- 
fer further discussion to a future publication. Parametric 
correlations at the mobility edge have not previously been 
investigated. We find if(t,r) ~ K exp(-ar|t| 1+? '/ d ), 
where a is a constant, giving 



R(0, A) - A T +^. 



(22) 



Thus the multifractal properties of critical eigenstates, 
via the exponent r], are reflected in parametric spectral 
statistics at the mobility edge. 



In summary, we have developed a simple and trans- 
parent approach for calculating spectral correlations in 
disordered metals. By averaging along a random walk 
through the ensemble of impurity configurations, we have 
expressed level statistics in all regimes in terms of the 
quantum return probability for a spreading wavepacket. 
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